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The MUNU experiment was carried out at the Bugey nuclear power reactor. The aim was the study of Vse~ 
elastic scattering at low energy. The recoil electrons were recorded in a gas time projection chamber, immersed 
in a tank filled with liquid scintillator serving as veto detector, suppressing in particular Compton electrons. The 
measured electron recoil spectrum is presented. Upper limits on the neutrino magnetic moment were derived and 
are discussed. 



1. Introduction 

The MUNU experiment was designed to study 
FeC" scattering at low energy, probing in partic- 
ular the magnetic moment of the neutrino. The 
detector was set up near a nuclear power reactor 
in Bugey (France) serving as antincutrino source. 
We report here on an analysis of data correspond- 
ing to 66.6 days of live time reactor-on and 16.7 
days reactor-off. 

2. Magnetic moments 

Many experiments now show that neutrinos 
have masses, and that the weak eigenstates Vi are 
superpositions of the mass eigenstates V£: 

yi = '^UtiVi i = 1,2,3; ^ = e,/x,T. (1) 

i 

Intense work is being done to determine more pre- 
cisely the mixings Uu and the masses mi. 

Besides masses, neutrinos can have magnetic 
moments. As shown in ref. |1I2I3| . the funda- 
mental magnetic moments are associated with the 
mass eigenstates, in the basis of which they are 



represented by a matrix /ijfe (j, fc = 1, 2, 3). Dirac 
neutrinos can have both diagonal and off-diagonal 
(transition) moments, while off diagonal moments 
only are possible with Majorana neutrinos. 

Astrophysical considerations put strong con- 
straints on Dirac neutrinos, and less strong ones 
on Majorana neutrinos (see j4l5j ). Upper limits 
applying in both cases from stellar cooling of or- 
der < 10~^^ have been derived. Magnetic 
moments large enough would lead to spin-flavor 
precession in the toroidal magnetic field in the in- 
terior of the sun, as discussed for instance in |6I7| . 
further complicating the oscillation pattern of so- 
lar neutrinos resulting from masses and mixings. 

These astrophysical bounds on magnetic mo- 
ments, however, depend on various assumptions 
and, to a certain extent, are model dependent. 
Direct measurements are under much better con- 
trol. So far the best limits stem from experiments 
studying V(,e~ or V^e" scattering, and looking for 
deviations of the measured cross-section from the 
one expected with weak interaction alone. The 
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cross-section is given by 
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with the contribution of the weak interaction in 
the first two Hnes, and that from the magnetic 
moment /Xg in the last one (Sj . Here E^ is the inci- 
dent neutrino energy, T the electron recoil energy, 
and the couplings are given by 
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The relative contribution of the magnetic moment 
term increases with decreasing neutrino and elec- 
tron energies. Therefore it is essential to have a 
low electron detection threshold, looking for neu- 
trinos from a low energy source. So far the sun 
and nuclear reactors have been used. Future ex- 
periments with radioactive sources are planned. 

The measured squared magnetic moment /i^ 
depends on the mixings and, in case of a large 
distance L from the source to the detector, on the 
propagation properties of neutrinos. For vacuum 
oscillations it is given by: 
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the momentum of Vk with mass ruk, for a given 
neutrino energy Ei,. In the case of matter en- 
hanced oscillations the dependence of the prop- 
agation eigenstates on the local density must be 
taken into account 

2.1. Solar neutrinos 

Here the source to detector distance is long, os- 
cillations are relevant and the measured quantity 



is denoted . Super-Kamiokande has measured 
with a very high statistical accuracy the electron 
recoil spectrum from the scattering of solar 
neutrinos ^ . The total rate is reduced because of 
neutrino oscillations, but the shape is seen to be 
in good agreement, within statistics, with that ex- 
pected assuming weak interaction alone. Beacom 
and Vogel looked to what extent excess counts 
at the low energy end from a magnetic moment 
can be ruled out. The assumption is made that 
such an excess is not compensated by a distortion 
due to oscillations. The limit < 1.5-10"^" fiB 
at 90 % confidence level (CL) was derived. 

2.2. Reactor neutrinos 

Nuclear reactors are strong sources of with 
energies ranging up to about 8 MeV. These are 
essentially produced in the beta decay of fission 
fragments, with however a significant contribu- 
tion at low energy from nuclei activated by neu- 
trons. Above 1.5-2.0 MeV the integral beta spec- 
trum of the fission fragments of the isotopes in 
which fission is predominantly induced by ther- 
mal neutrons (^^^U, ^^^Pu, ^'*^Pu) was mea- 
sured. These isotopes dominate. A simple proce- 
dure leads to the corresponding neutrino spectra, 
which are thus known well, with a precision of or- 
der 5 % or better ^Hl- In ^^^U, fission is induced 
by fast neutrons, and no such measurements are 
available. Less precise calculations (~ ±10 %) 
only can be used. But this isotope contributes 
only about 6-7 % in a conventional reactor. The 
sum spectrum above 2 MeV can be reconstructed 
with a precision of order ±5 %, knowing the rela- 
tive contributions of the various fissile isotopes at 



a given reactor. Studies of i^^p 



e^n scattering 



at reactors |10lll| . which sample the antineutrino 
spectrum above 2.5 MeV, show good agreement 
with these predictions. 

Below 1.5-2 MeV the neutrino spectrum can 
only be reconstructed from inherently less precise 
calculations as shown in ref . ■ In that work it 
is mentioned that neutron activations of the fissile 
isotopes, leading in particular to ^^^U, ^'^^Np and 
^^^U, and of the fission products, contribute. The 
yields have been estimated. 

In reactor experiments the source to detector 
distance is short compared to oscillation lengths. 
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and the magnetic moment searched for /i* is 
given by eq. Q with L set to zero. 

The Irvine group was the first to observe v^e" 
scattering 13 . The experiment was performed at 
the Savannah River reactor, with a 16 kg plastic 
scintillation counter surrounded by a Nal veto at 
11 m from the reactor core. A reactor-on minus 
reactor-ofF signal was seen, with a threshold at 1.5 
MeV. Analyzing the electron recoil data with the 
most recent knowledge of the Weinberg angle and 
the reactor spectrum, Vogel and Engel ,8, found 
a slight excess of events which can be explained 
by a magnetic moment /ig'"""* of order (2 — 4) • 

Also Vf,e~ scattering was studied at the Rovno 
reactor by a group from Saint Petersburg. The 
detector consisted of a stack of silicon sensors, 
with a total mass of 75 kg. A reactor-on minus 
reactor-ofF signal was seen above 600 keV, with 
a signal to background ratio of order 1:100. The 
limit < 1.9(1.5) • 10-1" at 95(68) % CL 

was reported p^j . 

More recently the TEXONO collaboration in- 
stalled an Ultra Low Background High Purity 
Germanium detector, with a fiducial mass of 1.06 
kg, near the Kuo-Sheng reactor in Taiwan jl5l . 
Here the approach is somewhat different. The 
threshold on the electron recoil (12 keV) is ex- 
tremely low. The reactor-on and reactor-off spec- 
tra were found to be identical within statistical 
errors. From that the hmit ^xf"''* < 1.3(1.0) • 
10-1° MB at 90(68) % CL was derived. 

3. MUNU 

In the aforementioned experiments the energy 
of the recoil electron candidate only was mea- 
sured. The MUNU collaboration |5ITB| has built a 
detector of a different kind, in which the topology 
of events is recorded. This allows a better event 
selection, leading to a lower background. More- 
over, in addition to the energy, the initial direc- 
tion of an electron track can be measured. A sec- 
ond parameter, the electron scattering angle, can 
therefore be reconstructed. This allows to look 
for a reactor signal by comparing forward elec- 
trons, having as reference the reactor to detector 
axis, with the backward ones. The background 



is measured on-line, which eliminates problems 
from detector instabilities, as well as from a pos- 
sible time dependence of the background itself. 




Figure 1. The MUNU detector at the Bugey re- 
actor. 



3.1. The detector 

The detector is described in details in ref . ^B] , 
and we only present here the essential features. 
The detector, made from radiochemically clean 
materials, was installed at 18 m from the core 
of a commercial reactor in Bugey (France) with a 
power of 2750 MWth. It emits neutrinos from the 
fission fragments of 235]j (•54 % qj-^ average over 
an annual reactor cycle), ^'^^Pu (33 %), ^'^^Pn (6 
%), and 238U (7 %). 

The central component of the MUNU detec- 
tor consists of a time projection chamber (TPC) 
filled with 3 bar of CF4 gas (figure P), acting as 
target and detector medium for the recoil elec- 
tron. CF4 was chosen because of its relatively low 
Z, leading to reduced multiple scattering, its ab- 
sence of protons, which eliminates backgrounds 
from Vf^p — > e^n scattering, its good drifting 
properties and its high density. As shown in fig- 
ure ^ the gas is contained in a cylindrical acrylic 
vessel of 1 m^ volume (90 cm in diameter, 162 
cm long, total CF4 mass 11.4 kg), which is en- 
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tirely active. The drift field parallel to the TPC 
axis (z-axis) is defined by a cathode on one side, 
a grid on the other side, and field shaping rings 
at successive potentials outside the acrylic vessel. 
The voltages were set to provide a homogeneous 
field leading to a drift velocity of 2.14 cm-/xs~^. 
An anode plane with 20 fim wires and a pitch of 
4.95 mm, separated by 100 /im potential wires, 
is placed behind the grid, to amplify the ioniza- 
tion charge. The integrated anode signal gives 
the total energy deposit. 

A pick-up plane with perpendicular x and y 
strips (pitch 3.5 mm) behind the anode pro- 
vides the spatial information in the x-y plane 
perpendicular to the z-axis. The spatial infor- 
mation along the z-axis is obtained from the 
time evolution of the signal. To reduce sys- 
tematics when comparing forward and backward 
events, the TPC was positioned orthogonally to 
the reactor-detector axis, which moreover coin- 
cides with the bisecting line between x and y 
strips on the pick-up plane. The anode wires are 
rotated by 45° with respect to the x-y plane. The 
TPC is thus absolutely symmetric between back- 
ward and forward directions with regards to the 
reactor-detector axis. 

The imaging capability of the TPC is illus- 
trated in figure 13 which shows an electron track. 
The energy resolution and calibration is obtained 
by comparing spectra measured using various 7 
sources (^^^Cs (662 keV), '^''Mn (835 keV) and 
22Na (1274 keV)) with simulations. The rela- 
tive energy resolution is found to be 8 % (1 cr) 
at 1 MeV, correcting for small variations across 
the anode plane. It scales with the power 0.7 of 
the energy, rather than with the square root. We 
think that this is because the electron attache- 
ment in CF4 in the strong field around the anode 
wires is rather high, affecting the statistics in the 
avalanche |16ll7j . The gain stability is monitored 
regularly throughout the data taking period, with 
sources and by measuring the spectrum of cosmic 
muons crossing the TPC. The data are corrected 
for small instabilities. The angular resolution is 
around 10° (1 a) at 1 MeV, for tracks scanned vi- 
sually, as derived from Monte-Carlo simulations, 
with a slight angular dependence jl8| . The energy 
dependence is only weak above 700 keV. 
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Figure 2. A 1 MeV electron in the TPC. The x- 
z and y-z projections are shown. The increased 
ionization which allows to determine the end of a 
track is clearly visible. 

The acrylic vessel is immersed in a steel tank 
filled with liquid scintillator and viewed by pho- 
tomultipliers, 24 on each side. It acts as a 
veto counter against the cosmics and as an anti- 
Compton detector to reduce the background from 
7's. The scintillator also sees the light produced 
by the avalanche around the anode wire, provid- 
ing an additional measurement of the ionisation 
charge. The primary light of heavily ionizing par- 
ticles such as a's is seen as well. A few a's are 
observed from remaining surface contaminations 
on the cathode. The primary light of minimum 
ionizing particles confined in the gas volume how- 
ever is below detection threshold. 

The TPC threshold is set at 300 keV, and that 
in the scintillator at 100 keV, with a minimum of 
5 photomultipliers hit. We note that the singles 
rates in the photomultipliers show no azimuthal 
dependence, which could result from hot spots. 

In normal data taking TPC events above 
threshold are read out, provided they are not in 
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coincidence with a 200 /iS signal started by scintil- 
lator pulses above 22 MeV. This eliminates direct 
cosmic hits, or neutrons associated with them. 
The 48 photomultipliers are read out first. The 
light from the photomultipliers on the anode side 
and the cathode side is compared. Reading pro- 
ceeds only if the relative difference is less than 
±30 %. Events corresponding to real tracks inside 
the gas volume were found to always fall within 
these limits . Discharges however can result in 
a larger imbalance, in which case they are read- 
ily eliminated. This helps reducing the average 
readout time. 



has the advantage of cleaner event selection. Also 
the risk of misidentifying the beginning of a track 
is much reduced. This approach is time consum- 
ing, however. To minimize the work load, that 
analysis was restricted to energies above 700 keV. 
Results from an improved version of the auto- 
matic scanning procedure will be discussed in a 
subsequent section. The data sets corresponds to 
66.6 days of live time reactor-on and 16.7 days 
reactor-off. 



4. Event selection 



Good events are single electrons contained in 
the TPC volume. The selection of neutrino scat- 
tering events proceeds in two steps. First an au- 
tomatic filtering eliminates obviously bad events, 
namely events 

• identified as a's or discharges from their 
topology (high ionization in a small vol- 
ume), 

• or in delayed coincidence (80 /iS, corre- 
sponding to the TPC length) with a signal 
in the scintillator, 

• or events not contained in a fiducial volume 
of 42 cm radius, 

• and finally events with a fast rise time, due 
to particles crossing the amplification gap 
between the pick up plane plane and the 
grid, in either direction. Here the avalanche 
light signal is used, as described in ref. ^Hl • 

The precise live time is derived in this process on 
a daily basis. It is found to fluctuate around 65 
%. It is limited primarily by the total veto time 
of the scintillator (11 %) and the dead time of 
the TPC itself, caused by the relatively long data 
read-out and data transfer time (24 %). 

Then a final scan is performed. In a separate 
publication, we reported first results from a crude 
automatic procedure |2()j , carried out with a pat- 
tern recognition program. Here we present an 
analysis based on a visual scan of events, which 
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Figure 3. Distribution of Oreac for single electron 
events, top reactor-on, bottom reactor-off, visual 
scanning. 
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In this second scan both the x — z and y — z 
projections of an event are scrutinized, as well as 
the evolution in time of the anode and scintilla- 
tor signals. Continuous electron tracks only are 
retained. The end of the track is identified from 
the increased energy deposition, due to the higher 
stopping power, as exemplified in figure[51 Events 
with a second high deposition along the track, in 
particular near the other end, are discarded. This 
improves the background suppression, at the cost 
of events with a delta electron, which however 
contribute negligibly. 

Then the tangent at the start of the track is 
determined by eye. From that the angles Oreac 
with respect to the reactor core-detector axis is 
determined, as well as the angle Odet with respect 
to the TPC axis. 

An excess of events from the anode side was 
observed. It is presumably due to additional ac- 
tivities, resulting from the greater complexity of 
the readout system, and to a larger inactive vol- 
ume in the scintillator because of the stronger and 
thicker acrylic lid. For that reason only electrons 
emitted in the half sphere from the cathode side 
{ddet < 90") are accepted. This reduces the ac- 
ceptance by a factor 2, but leads to a better signal 
to background ratio. 

5. Results from the visual scan 

In figure |31 we show the distribution of 
cos{9reac) of single contained electrons, for both 
reactor-on and reactor-off. The slightly non linear 
angular response of the TPC, and its geometry, 
explain the accumulation of events at cos{6reac) 
around 1,-1 and for reactor-off. The distribu- 
tion is however identical in forward (cos (^reac) — 
1) and backward {cos{0reac) — ^1) directions. 
The reactor-on spectrum shows a clear excess of 
events in forward direction from T'ee~ scattering. 
Uncertainties from instabilities of the gain, the 
veto rate, the live time, cancel out in this forward 
minus backward comparison. 

For each electron event the neutrino energy E^, 
is reconstructed from the electron recoil energy T 
and from the scattering angle, taken as 9reac- For- 
ward event candidates are defined as those with 
positive neutrino energy: E^, > 0. To select back- 



ward events the neutrino energy which must be 
positive is that calculated using tt — fl^eac a-s scat- 
tering angle. This procedure has an acceptance 
close to 100 % above 700 keV, as determined by 
Monte-Carlo simulations. It takes into account 
that electrons are emitted in an narrower cone 
whith increasing energies, and is somewhat more 
sophisticated than the application of a crude cut 
on cos{0reac) ■ In practicc, however, with a thresh- 
old of 700 keV, it produces almost the same re- 
sult as the application of a cut cos{0reac) > 0.7 

or C0s{9reac) < "0.7. 

The energy distributions of both forward and 
backward events are displayed in figure^ A clear 
excess of forward events (458 in total) over back- 
ward events (340) is seen. The total forward mi- 
nus backward count rate above 700 keV is thus 
1.77±0.42 day~^. The background, given by the 
backward events, shows a steep low energy com- 
ponent ending at about 1.2 MeV. Above it is seen 
to be fairly low, and much flatter. The back- 
ground was observed to increase slightly during 
the course of the experiment, by about 10 % in 6 
months, possibly because of outgassing. A small 
drop in efficiency of the anti-Compton, even with- 
out an observable reduction in count rate, could 
also explain this. 

As a cross check the reactor-off recoil spec- 
tra, both forward and backward, measured af- 
ter the reactor-on period, were reconstructed as 
well, and are displayed in figure [S] They are in- 
deed identical within statistics, the integrated for- 
ward minus backward rate above 700 keV being 
-1.02 ± 1.00 day-^ 

The difference of the forward minus backward 
reactor-on spectra is shown in figure IHI The ex- 
pected event rate was calculated using the best 
knowledge of the reactor spectrum as described 
above, and taking into account the known activa- 
tions of fissile isotopes and fission products. The 
uncertainty is around 5 % above 900 keV, as dis- 
cussed previously, and larger below. The various 
acceptances of the event selection procedure were 
determined from measurements with sources, and 
by Monte Carlo simulations using the GEANT3 
code. Data sets including recoil electrons from 
neutrino scattering were produced, filtered and 
scanned just as the real data. The containment 
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Figure 4. Energy distribution of forward and 
backward events, reactor-on, visual scanning. 



16.7 days 
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Figure 5. Energy distribution of forward and 
backward events, reactor-off, visual scanning. 



efficiency in the 42 cm fiducial radius for recoil 
electrons was found to vary from 63 % at 700 keV, 
50 % at 1 MeV to 12 % at 2 MeV. The relative 
uncertainty is of order 2-3 %. Some tracks with 
weird topologies cannot be reconstructed, reduc- 
ing the acceptance by 4 %, relatively speaking. 
Similarly the remaining acceptances together, in- 
cluding that of the > cut, lead to an ad- 
ditional 6 % reduction. Neutrino interactions in 
the copper cathode give rise to electrons which 
can escape into the gas volume. These cannot 
be vetoed. Taking into account all cuts they in- 
crease the expected rate between 700 keV to 2 
MeV by some 3 %. The relative uncertainty on 
the global acceptance is of order 7 %, leading to 
a total uncertainty of 9 % on the total expected 
rate. 

The total expected rate above 700 keV as- 
suming a vanishing magnetic moment was found 
to be 1.02±0.1 day~^, to be compared with 
1.77±0.42 day~^ measured, as mentioned above. 
There is thus a certain excess of measured events, 
which has however only a small statistical sig- 
nificance. The calculated energy distribution is 
shown in figure The excess counts are seen to 
be in the region below 900 keV. 

To be more quantitative, a "^^^ calculated 
using 100 keV bins as shown in figure from 700 
keV to 1400 keV, and then a bin from 1400 to 




Figure 6. Energy distribution, forward minus 
backward events, reactor-on, visual scanning. 



2000 keV. Gaussian statistics applies then to all 
bins. The error on the expected rate is small in 
comparison to the statistical uncertainties, and 
turns out to be neghgible. The was calculated 
this way for several values of the squared mag- 
netic moment (/ig'*"''*)^, constraining it to the re- 
gion = (0.77 ± 0.92) ■ 10-20 ^1%. This 
is consistent with a vanishing magnetic moment. 
Renormalizing to the physical region {nf^"^*)"^ — 



we find the limit < 1.4(1.2) • 10"!° fis at 

90(68)% CL. 

Nevertheless even the best fit is not very sat- 
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isfactory. The (9-9 for 7 degrees of freedom) 
is on the high side. More troublesome, as seen in 
figure ini the inclusion of a magnetic moment in 
the calculated spectrum improves the agreement 
with the data in the two first bins from 700 to 
900 keV, but makes it worse in the upper ones. 

These two first bins are solely responsible for 
the high total rate. We note that such an ex- 
cess is also visible, unfortunately also with limited 
statistical precision, and appearing at somewhat 
smaller energies, in the data of ref. This ex- 
cess may well result not from a magnetic moment, 
but instead from sources not taken into account 
when evaluating the reactor neutrino spectrum. 
As mentioned above, the low energy part, which 
contributes significantly to the electron spectrum 
below 0.9 MeV, is not so well known. Additional 
neutron activations may contribute, beyond the 
known ones of the fissile isotopes and the fis- 
sion products. In that sense it seems safe at 
this stage to restrict the analysis to electron en- 
ergies above 900 keV. There the measured event 
rate is 0.41 ± 0.26 day^^, in good agreement with 
the expected one 0.62 ± 0.05 day^^. From the 
event rate and the energy distribution, using the 
same prescription as above, the allowed range 
(^^shoTt^2 ^ (_o.95 ± 0.95) • 10-20 ^1% is found. 
The best is 1.51 for 5 degrees of freedom. This 
yields the limit 

^short ^ 10(0.8) • 10-1° at 90(68) %CL (4) 
somewhat more stringent than the previous one. 

6. Automatic scanning 

Extending the visual scan to energies below 
700 keV is not feasible because of the rapid in- 
crease of the background below that energy, lead- 
ing to an unmanageable workload. But the auto- 
matic scanning procedure mentioned in j20| has 
been improved, and used on the same data set 
|19| . The program first identifies electron tracks, 
searches the vertex, and then fits the beginning 
of the track. 

Electron tracks as identical as possible to real 
tracks were produced by Monte-Carlo. The elec- 
tronics noise was included. The tracks were an- 
alyzed with the same program to determine the 



global acceptance, as well as the angular resolu- 
tion. The angular resolution is inferior to that 
obtained in the visual scanning, varying from 31" 
at 300 keV, to 18° at 700 keV and finally 12° 
at 2 MeV. This reduces the acceptance of the 
El, > cut in particular at low energy. This, 
combined with poorer electron track identifica- 
tion, leads to an acceptance of 20 % at 300 keV, 
increasing to 54 % at 1 MeV and remaing con- 
stant above that energy. Events with 9det < 100° 
were taken. The other cuts and acceptances are 
the same as in the visual scanning. The back- 
ground suppression is inferior by a factor 3 or 
more, so that the statistical precision is less. But 
the method can be applied to energies extend- 
ing down to 300 keV. We show in figure [7| the 
forward minus backward spectra for reactor-on. 
It is in good general agreement with the spec- 
tra from the eye scan. The high bin at 750 keV 
is reproduced. But, admittedly with relatively 
large statistical errors, a catastrophic rise when 
going to lower energies can be ruled out. Apply- 
ing the same procedure as above to the data one 
obtains the limit /if < 1.7 • 10^1° ^J.B at 90 
% CL. It does not change with the threshold in 
the range 300 to 800 keV, the better sensitivity 
at low energy being offset by the larger statistical 
uncertainties. 
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Figure 7. Energy distribution, forward minus 
backward events, reactor-on, automatic scanning. 
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7. Conclusions 

The MUNU experiment studied Vf,e~ scatter- 
ing at low energy near a nuclear reactor. For elec- 
tron recoil energies above 900 keV, good agree- 
ment is seen with expectations assuming weak 
interaction alone. From this the following limit 
on the magnetic moment of the neutrino can be 
derived: /^f < 1.0(0.8) • 10-^° /xb at 90(68) % 
CL, limited primarily by statistics. 

In any event this bound rules out the possible 
indication for a relatively large magnetic moment 
from the Savannah River experiment |S] , and im- 
proves on the limit from ref. |14j and |15j . It is 
somewhat more stringent than the limit on /i*"' 
from solar neutrinos 3 , which however does not 
apply to exactly the same quantity. 

The authors want to thank Edf-CNPE Bugey 
for the hospitality. This work was supported by 
IN2P3, INFN, and SNF. 
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